Quantum control and entanglement in a chemical compass 
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The radical pair mechanism is one of the two main hypotheses to explain the navigability of 
animals in weak magnetic fields, enabling e.g. birds to see the Earth's magnetic field. It also plays an 
essential role in the field of spin chemistry. Here, we show how quantum control can be used to either 
enhance or reduce the performance of such a chemical compass, providing a new route to further 
study the radical pair mechanism and its applications. We study the role of quantum entanglement 
in this mechanism, and demonstrate intriguing connections between radical-pair entanglement and 
the magnetic field sensitivity of the compass. Beyond their immediate application to the radical 
pair mechanism, these results also demonstrate how state-of-the-art quantum technologies could 
potentially be used to probe and control biological functions. 



Introduction. — It is known that many species, includ- 
ing birds, insects and mammals, use the Earth's mag- 
netic field for orientation and navigation[l|. To explain 
this remarkable ability, two main hypotheses have been 
proposed: a magnetite-based mechanism and a radical 
pair biochemical reaction mechanismjlj. Since the radi- 
cal pair mechanism (RPM) was first proposed in pioneer- 
ing work by Schulten et al. [1] , a chemical compass model 
for migratory birds, based on such a mechanism 0] has 
been widely studied. Evidence suggests that the RPM 
is indeed linked to the avian magnetoreception 0, • It 
was recently demonstrated in spin chemistry experiments 
that a photochemical reaction can act as a compass even 
in a magnetic field as weak as the geomagnetic field Q. 
The underlying mechanism in such a chemical compass 
is clearly of quantum mechanical nature. However, the 
detailed role of quantum interactions, in giving rise to 
entanglement and (de-) coherence, are little understood 
0. On the other hand, one can observe growing interest 
in the role of quantum coherence for biological processes 
in general [§], and specifically in photosynthesis [§]. A 
deeper understanding of the role of quantum mechan- 
ics in biology will eventually come along with the ability 
to control biological processes at the level of individual 
molecules. In physics, various kinds of quantum control 
techniques have been developed, specifically in the field 
of quantum information processing and quantum metrol- 
ogy 13, III- The question thus naturally arises to what 
extent these or similar techniques could be applied to 
test and refine certain biophysical hypotheses, such as 
the chemical compass model for animal magnetorecep- 
tion? Can we use quantum technologies that have pri- 
marily been developed to control man-made microscopic 
systems, to study the behavior of living things — e.g. 
birds, fruit files, or plants — in a detectable way? 

In our work, aiming at the above questions, we will 
revisit the RPM and the chemical compass model using 
concepts and techniques from quantum information. The 
RPM can serve both as a magnetometer or as a compass, 
depending on the molecular realization. For simplicity, 
we will refer to both cases as "compass" in the following. 



First, we demonstrate that quantum control ideas can be 
applied to experiments in spin chemistry and potentially 
also to study the magnetoreception of certain animals. 
We propose several quantum control protocols that can 
be used to either enhance or suppress the function of a 
chemical compass. Assuming that the model provides 
the correct explanation for magnetoreception of certain 
species, we predict that they would loose or regain their 
orientability in appropriately designed experiments using 
such quantum control protocols - given that such experi- 
ments could be carried out safely. Our calculations show 
that the RPM can not only detect weak magnetic fields, 
but it is also sensitive to quantum control even without 
the presence of a static magnetic field. These results offer 
a new means to study experimentally the RPM, also in 
comparison with other conceivable mechanisms such as 
those in man-made magnetometers Bill. 

Second, we investigate whether entanglement is a nec- 
essary ingredient in animal magneto-reception, which 
seems appealing in the light of the important role this 
concept has gained in fundamental discussions on quan- 
tum mechanics and its wider implications. As the sensi- 
tivity of the chemical compass depends on the initial state 
of the radical pair, it is natural to ask whether it needs to 
be quantum mechanically entangled - thereby excluding 
any conceivable classical mechanism - or whether clas- 
sical correlations would be sufficient. We find that the 
answer largely depends on the radical pair lifetime. For 
specific realizations of the RPM, e.g. those in recent 
spin-chemistry experiments [l^ . entanglement features 
prominently and can even serve as a signature of the un- 
derlying spin dynamics. However, when the radical pair 
lifetime is extremely long, as it is believed to be the case 
in the molecular candidate for magneto reception in Eu- 
ropean robins [isl [l6j . entanglement does not seem to 
play a significant role. 

Radical Pair Mechanism. — We consider a photochem- 
ical reaction that starts from the light activation of a 
photoreceptor, followed by an electron transfer process; 
two unpaired electrons in a spin-correlated electronic sin- 
glet state are then carried by a radical pair. The effective 
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environment of a radical pair mainly consists of their indi- 
vidual surrounding nuclei. The Hamiltonian of a radical 
pair is of the form [17 



H 



J2 Hk = -JeB-J2Sk 



fe=1.2 



where 7e = —gefJ-B is the electron gyromagnetic ratio, 
Xkj denote the hyperfine coupling tensors and Sk, Ikj 
are the electron and nuclear spin operators respectively. 

The initial state of a radical pair is assumed to be the 
singlet state |§) = -^(Iti) — lit))i which subsequently 
suffers from de-coherence through the hyperfine interac- 
tions with the environmental nuclear spins. The initial 
state of the nuclear spins at room temperature can be ap- 
proximated as pfc(O) = ^j^j/dj, where dj is the dimen- 
sion of the jth nuclear spin. The charge recombination 
of the radical pair goes through different channels, de- 
pending on the electron-spin state (singlet or triplet). In 
particular, the yield of products formed by the reaction 
of singlet radical pairs can be calculated as [l3| 



<^s{t) = / r,{t)f{t)dt 



(2) 



where rdt) is the radical re-encounter probability dis- 
tribution, and f{t) ~ {§\ps{t)\§) is the fidelity between 
the electron spin state Ps{t) at time t and the singlet 
state. The ultimate activation yield $s = ^s{t ^ oo) in 
cryptochrome is believed to affect the visual function of 
animals 

We have followed the established theory for the dy- 
namics of the RPM IT], |l8| and computed the full quan- 
tum dynamics of the combined system of electron spins 
and nuclear spins. Technically, we employ the Chebyshev 
polynomial expansion method [lit to numerically calcu- 
late the exact evolution operator Uk{t) = exp{—iHkt), 
and thereby all relevant physical quantities. We first con- 
sider the well-studied photochemical reaction of pyrene 



(Py-/iio) and N,N-dimethylaniline (DMA-Zin) UJ], for 
which the hyperfine couplings are isotropic |14{ . and the 
tensor A^^ in ([1]) simplifies to a number X^^ ■ We study 
the role of entanglement in this radical pair reaction and 
propose new experiments based on quantum control. We 
then generalize our results to the cryptochrome radi- 
cal pair of FADH'-O' , which is the molecular candi- 
date believed to be involved in avian magnetoreccption 
[isl [Tif . We thereby show that our protocols work also 
for anisotropic hyperfine interactions, which are essential 
for direction sensitivity of the magnetic field ^ . 

Magnetic Field Sensitivity under Quantum Control. — 
The magnetic-field sensitivity A of the radical pair re- 
action [Py-h'{^ DMA-/ij"|] is quantified by the derivative 
of the activation yield with respect to the magnetic field 
strength B [l3], i.e. 



MB) 



5$. 



(3) 



We assume that the external magnetic field points in the 
z direction. The key ingredient in the RPM are the hy- 
perfine interactions, which induce a singlet-triplet inter- 
conversion (mixing) depending on the magnetic field [l7j . 

as an example 



Using an exponential model r^t) 



-kt 



for the re-encounter probability distribution [17| , we plot 
in Fig. [ija) the magnetic-field sensitivity A as a func- 
tion of B. Our numerical simulation agrees well with the 
experimental results in 
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Figure 1: Magnetic field sensitivity A of the radical pair reac- 
tion [Py-/ijQ DMA-/i^"5^] as a function of the magnetic field 
B. (a) N: Singlet initial state; Z: under Z control; KB 
(RB-X): alternating magnetic field without (with) X con- 
trol, (b) N: Singlet initial state; Tq: Triplet initial state 
|To); Sep: Optimal sensitivity for separable initial states; 
CS-P: applying a -l-X pulse on the initial separable state 
pc = (Iti) (tJ-l + 14-t) (4-tj)/2- The recombination rate con- 
stant is fc = 5.8 X 10* s~^ [Ij]) and the control time is Tc — 0.5 



Studying the performance of the radical-pair mecha- 
nism under quantum control would allow us to test the 
role of entanglement and further details of the RPM 
in spin chemistry experiments. As a simple example, 
consider a periodical pulse sequence with 7r-pulses ap- 
plied at times t = niTc along the z direction; the effec- 
tive Hamiltonian to the first order is given by = 



Such kind of con- 



trol can actually enhance the performance of quantum- 



coherence based magnetometers, see e.g. [ij, ll3[- How- 
ever, in case of the RPM, the magnetic-field sensitivity 
becomes greatly suppressed, as can be seen in Fig.[TJa). 
We can show that, whenever one applies more general 
decoupling protocols to promote quantum coherence in 
a radical pair reaction, its magnetic-field sensitivity will 
generally be reduced [20|]. This demonstrates that it is 
in fact the decay of coherence, i.e. rfe-coherence, rather 
than coherence itself, that plays an essential role for the 
magnetic-field detection in RPM, different from the situa- 
tion in magnetometers using e.g. NV-centers in diamond 
11, [11. 

To demonstrate a potentially positive effect of quan- 
tum control on a chemical compass, we consider a situa- 
tion where the magnetic field alternates its direction pe- 
riodically at times t = ruTa which will disturb the proper 
functioning of compass. (This situation is reminiscent of 
an experiment with birds in an oscillating field [^, even 
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though the cause of the compass disfunction is here dif- 
ferent.) If we now apply tt-X pulses at the same times 
t = iTiTa, the chemical compass will recover its function 
as the transitions between |S) and |T±) induced by the 
residual xx hyperfine interactions are still affected by the 
magnetic field [see Fig. [TJa)]. 

Entanglement and Magnetic Field Sensitivity. — We 
have hitherto assumed, as is usually done, that the radi- 
cal pair starts in a perfect singlet state, i.e. that quantum 
coherence is fully maintained during the pair creation. In 
reality, the initial state Ps{Q) of radical pairs will never 
be a perfect singlet (i.e. pure state), but a mixed state 
with a certain singlet fidelity /(O) = (§1/3,5(0) |§) < 1. It is 
known that the value of /(O) has to be sufficiently close to 
unity, otherwise the state may also be described by clas- 
sical correlations. We therefore ask: Is entanglement, as 
a genuine quantum signature, needed at all to account 
for the efficiency of the magnetic compass? Or could the 
latter be explained by mere classical correlations? To an- 
swer this question, we have randomly chosen 5000 differ- 
ent initial states from the set of separable states and cal- 
culated the maximal achievable magnetic field sensitivity 
for every value of B, see Fig.[IJb) (o). We find that, in the 
operating region of the compass (around B = 4mT), the 
maximal achievable sensitivity for separable states stays 
significantly below the sensitivity for the singlet state, 
and this maximum sensitivity is in fact attained by the 
classical mixture pc = (|ti) (til + lit) (itl)/2. In turn, if 
nature is allowed to optimize the initial state from the full 
set of states, including the entangled states (e.g. |§) and 
|To), see Fig.[ljb)), the optimum magnetic-field sensitiv- 
ity will typically be much higher than for any separable 
state. On these grounds one can say that entanglement is 
indeed helpful, and it is specifically entanglement rather 
than mere quantum coherence. 

To test experimentally whether the initial state of 
the radical pair is indeed (close to) a singlet state, i.e. 
whether de-coherence can be neglected during the radi- 
cal pair creation, one could apply a ^-X pulse as the re- 
action starts. For an initial singlet state, which remains 
invariant under such pulse, the magnetic field sensitivity 
will remain unchanged, whereas for an initial classical 
mixture it will collapse, see Fig. [IJb). 

As entanglement seemingly plays a role in the RPM 
with Py-DMA, we have studied its dynamics and its 
quantitative connection to the magnetic-field sensitiv- 
ity. Similar to the activation yield, we define = 
rc{t)E(t)dt to quantify the effective amount of entan- 
glement that is present in the active radical pairs during 
the reaction, where E{t) is chosen to be the entanglement 
measure of concurrence [2l| at time t. The first derivative 
with respect to the magnetic field, A^ = 8^^/93, quan- 
tifies how sensitive this effective entanglement is with re- 
spect to variations of the magnetic field. 

In Fig. [2ja), we see that A^ and A are correlated in 
the regions of I and II, displaying monotonic relations 
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Figure 2: Connection between quantum entanglement and 
magnetic field sensitivity in the radical pair reaction [Py- h'^ 
r>MA-h'{^]. (a) Sensitivity of effective entanglement A_b vs. 
sensitivity of singlet yield A. The recombination rate con- 



stant is k = 5.8 X 10 s 



The blue arrows indicate 



variation of A_e and A when the magnetic field changes from 
B = 0.5 mT to B = 8 mT. (b) Discontinuity of the lifetime 
of entanglement Te as a function of B. 



with different linear ratios. This result is remarkable in- 
sofar as that the time during which entanglement exists 
is significantly shorter than the reaction time T^. for the 
value oi A{B,t) d^{t)/dB to saturate However, 
it can also be seen from Fig. [5{a) that A e changes dra- 
matically at the crossover between regions I and II. This 
step-like behavior relates to the discontinuity of the en- 
tanglement lifetime T^ = max{t|£'(<) > 0} as the mag- 
netic field increases, see Fig. [UJb). In the region of I, 
Tb is much shorter than the reaction time T^, while it 
jumps to a larger value comparable with T,. during the 
crossover from the region of I to II. When we further 
increase the magnetic field, T^; exhibits more kinks but 
with less increment. This effect originates from the finite 
size of the nuclear spin bath of the electron spins, and 
is a clear signature of the system-environment dynamics 
underlying the RPM 

Applications to Animal Magnetoreception. — In order 
to account for a direction sensitivity of the singlet yield, 
which is necessary for compass function, the hyperfine 
couplings must be anisotropic 0|. Here we consider an 
example of such a radical pair, FADH'-O*", which was 
proposed as a likely molecular candidate underlying the 
magnetoreception of European robins [H, [l3| , but it may 
also play this role in other species. 

The direction of the magnetic field in ([T]) with respect 
to the reference frame of the immobilized radical pair is 
described by two angles (0,0), i.e. B = B{sm9 coscf), 
sin 6* sin (/), cos 6*) . Without loss of the essential physics, 
we here assume that = 0, and investigate the depen- 
dence of the singlet yield $s on the angle 9 when we 
apply quantum control. First, it can be seen from Fig. [3] 
(a) that the angular dependence of the singlet yield is 
much suppressed if one applies 7r-pulses along the same 
direction as the magnetic field, which can distinguish the 
RPM from other potential mechanisms for magnetore- 
ception [13, 13 1 . Next, we study the scenario that the 



magnetic field changes it direction periodically at times 
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Figure 3: Singlet yield $s of the radical pair reaction [FADH*- 
O*"] as a function of the angle 6 with the magnetic field 
S = 46 ^T. (a) N: Without quantum control; Qc: Apply- 
ing TT-pulses along the direction of the magnetic field, (b) 
RB (RB-X): Effect of an alternating magnetic field, without 
(with) additional quantum control pulses perpendicular to the 
direction of the magnetic field. For comparison, the RB curve 
has been shifted downwards by 0.1. The recombination rate 
constant is fc = 5 x 10^ and the control times are Tc ~ 10ns 
and Ta ~ 10ns. 



t = riTa as in the previous section. As expected, the 
angular dependence is again greatly suppressed, as can 
be seen from the lower curve in Fig. EJb). However, if 
one applies 7r-pulses perpendicular to the direction of the 
magnetic field this will re-induce an angular dependence, 
see Fig. |3l^b). Furthermore, we find that even without a 
static magnetic field, quantum control can induce an an- 
gular dependence of the singlet yield as shown in Fig. 2] 
In other words, if one would be able to design a behav- 
ior experiment with animals that use a chemical compass 
to sense the magnetic field, in such an environment they 
would lose or regain their orientability, depending on the 
applied control fields. This would provide further evi- 
dence for the RPM as the underlying mechanism, and it 
could help to narrow down the possible candidates of rad- 
ical pairs in animal magnetoreception. It is however not 
clear how such experiments could be carried out safely. 
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Figure 4: Singlet yield $3 of the radical pair reaction [FADH*- 
0*~] as a function of the angle 6, relative to the direction in 
which the vr-pulses are applied. There is no external static 
magnetic field, i.e. B = Q. The recombination rate constant 



is fc = 5 X 10 s , and the control time is Tn 



100 ns. 



Different from the example of Py-DMA we find that 
here the entanglement only exists in a time range 
10 ns) that is much shorter than the expected radical 



pair lifetime ('-^ 2 — 10/is). We would thus not expect 
entanglement to play a significant role in this context. 
To check this further, we have computed the achievable 
sensitivity of the compass for different initial states and 
found that a substantial part of all separable states can 
account for an angular dependence that is as high as 
(or even higher) than for the singlet state [2^. This 
means that - in contrast to Py-DMA - the radical-pair 
entanglement does not seem to be a necessary ingredient 
for a chemical compass based on FADH'-O'". 

Summary and Outlook. — We have demonstrated how 
quantum control can influence radical pair reactions and 
the function of a chemical compass. The presented proto- 
cols can in principle be applied to existing spin chemistry 
experiments, even though the implementation of coherent 
spin control [23, 24 1 in this context needs to be further 
developed. They might also provide a route for future 
experiments with biological systems (including animals 
or plants) that are expected to exploit the RPM; in this 
case a much more careful study would be required, in par- 
ticular regards the potential side effects of short control 
pulses on biological tissue. 

We found interesting connections between entangle- 
ment and the magnetic field sensitivity when the radical 
pair lifetime is not too long compared to the coherence 
time. Otherwise, the role of coherence and entanglement 
seem to be insignificant. Whether or not birds or other 
animals use entanglement for their ability to orient them- 
selves in the earth magnetic field remains an open ques- 
tion, whose answer will depend on the specific molecular 
realization of their chemical compass. 

As a bio-mimetic application of practical relevance, it 
would be interesting to explore the possibility of simulat- 
ing a radical-pair mechanism in more controllable quan- 



tum systems, such as NV centers in diamond |12l. 1131.12 



to design an ultra-high fidelity sensor for the detection of 
weak fields or forces. 
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Supplementary Material 

This is supporting material for our paper. We derive the completely positive map for the dynamics of a central spin 
coupled to its surrounding nuclear spins with isotropic hyperfine interactions, and use it to investigate the evolution 
of entanglement in the radical pair reaction. To identify the role of entanglement in the radical pair mechanism, we 
randomly choose the initial state from the set of separable states, to compute the optimum magnetic field sensitivity 
for these states and compare it with the sensitivity for an inital singlet state. Further details arc provided to clarify 
the connections between quantum entanglement and magnetic field sensitivity. To illustrate the essential role of 
the (de-phasing) nuclear spin environment in a chemical compass, we investigate a hypothetical reference model of 
bosonic thermal bath and compare it with the present results. 

Molecular structures of radical pairs. — The molecular structures for the radical pair Py-DMA are displayed in Fig.[TJ 
Py-^io has ten spin-^ hydrogen nuclei, while the DMA-Zin has eleven spin-^ hydrogen nuclei and one spin-1 nitrogen 
nucleus; the nuclear spin of carbon is 0. In our numerical simulations, without loss of essential features, we have 
considered three groups of equivalent nuclei in each radical that have the largest hyperfine couplings as in [l| , i.e. the 
radical Py-/iio interacts with ten spin-i surrounding nuclei with the hyperfine coupling constants A^iJ' = 0.481 mT 
(4x H), A^-^' = 0.212 mT (4xH), A^-J' = 0.103 mT (2xH) [l|, and the radical DMA-Zin is dominantly coupled with 

seven spin-i nuclei with A^ ^ = 1.180 mT (6xH), A^f = 0.520 mT (IxH), and one spin-1 nucleus with xfj = 1.100 
mT (Ix N), see Table 2 of 




Figure 1: (Color online) Molecular structures of the radical pyrene (Py-/iio) (left) and N,N-dimethylaniline (DMA-ftn) (right). 
Green: Carbon; Grey: Hydrogen; Blue: Nitrogen. 

The flavin radical FADH* is displayed in Fig. [51 We consider the dominant hyperfine couplings from two spin-1 
nitrogen nuclei and three spin-i hydrogen nuclei as in The superoxide radical 0|~ is devoid of the hyperfine 
couplings, which is likely to lead to higher sensitivity against weak magnetic fields 0, 




Figure 2: (Color online) Molecular structure of the flavin radical FADH*. Green: Carbon; Grey: Hydrogen; Blue: Nitrogen; 
Red: Oxygen. 

Completely positive map for electron spin dynamics. — Wc here derive the completely positive map for the case of 
isotropic hyperfine interaction. The Hamiltonian for a central unpaired electron spin coupled with a nuclear spin bath 
is written as 

H, = fUbS, + ^kS ■ I^^-) (4) 
fc 

where mt ~ —jeB. The presently available theories for the central spin problem usually resort to the perturbation 
approach, based on certain approximations, e.g. the quasi-static approximation or the limit of large magnetic fields 
and/or large spin bath polarizations. For our present purpose, these approximations arc only of limited use since, in 
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the radical pair mechanism, one is particularly interested in the region of low fields, and the number of most relevant 
surrounding nuclei is ~ 10, in contrast with ^ 10^ in quantum dots. 

It is straightforward to show that the total angular momentum of the electron and nuclear spins, Mz = Sz + Iz, 
where = J2k ^^''^ ' conserved for the Hamiltonian in Eq. (g]), i.e. [Mz,Hc] = 0. By introducing as the basis 

of eigenstates of /z, i.e. hW^) = '^l'*?™), where n labels the eigenvalues and m is a degeneracy index, we can express 
the initial state of the spin bath as P(,(0) = {^j, Ik/dk = ^ X^n m \v™){'^P^\^ where d = J^j, dk is the total dimension of 
all the (relevant) nuclear spins. Thus, under the coherent evolution Uc ~ exp{—itHc), the joint state of the central 
spin and the nuclear spins evolves as 

-> n)\^ln) + \i)\V>lnn) (5) 

^ It) b,7j.> + 1;> i^L) (6) 

wherein ||) and {l\ denote the eigenstates of Sz = ^cr^, and Iflnn) belongs to the eigenspace of Iz associated to the 
eigenvalue n + i. The fact that the total angular momentum is conserved results in orthogonality relations for the 
nuclear spin states: 

\^ij,\^ij±\^-^^)±\^ij (7) 

The inner products of these vectors are zero, as they belong to orthogonal subspaces (or are null vectors). By recalling 
the notation ^ = |t) (t|, ^ = \i) = It) (tl, = \i) (t|, we obtain 

= Tr^lO(^^J = 0(8) 

71, m 

in which we have used the relation in Eq. ([7|). In a similar way, one can show that /ig- = yUi± = /i±o = A'±i = M++ = 
/i = 0. Moreover, it is easy to verify that 

Moo = Mil = ^ + -^Tr {UcOzUlcTz) 

Thus, the dynamics of the central spin, which is calculated by tracing out its spin bath degrees of freedom as 
Ps{t) = Trt{e"*^'=*[ps(0) ® pt(0)]e*^'=*}, can be exphcitly expressed as 

m ■■ It) (tl ^ at It) (tl + (1 - at) \i) (;| 

i;)ai ^ (i-aoit)(ti+«ti;)ai 
It) ai ^ It) ai 
it)(ti ^ «nt)(ti 

from which we can obtain the completely positive map for the central spin dynamics. If we rewrite the evolu- 
tion operator Uc = cxp{—itHc) in the form of Uc = '^fj_ ^IfJ^cii^l ® Ufj_^, we get the above dynamic parameters 

at = Tr([/oo[/oo)' ^nd Kt = Tt{UqqUIj^). 

Dynamics of quantum entanglement in Py-DMA. — In Fig. [31 we plot the evolution of entanglement between two 
unpaired electron spins in the radical pair reaction [Py-h'^ DMA-h'-^]. We have used the concurrence Q as the 
measure of two-spin entanglement, which vanishes on separable states and assumes its maximum value, 1, on maximally 
entangled states such as the singlet state. We also compare the exact value of entanglement with the estimated best 
lower bound of entanglement e{t) = infp{£'(p)|Tr(p'Ps) = fs{t)} = max{0, 2fs{t) — l}. The agreement between them is 
good (even though not perfect). This fact supports our statement about how one could possibly estimate the amount 
of entanglement from experimentally accessible information. 

We also plot the dynamics of entanglement under X and Z control in Fig. [3] It can be seen that entanglement 
survives indeed for a longer time if quantum control is applied. 

Protecting coherence is not helpful. — In the quantum coherence based magnetometer, e.g. with NV centers in 
diamond [fl 0], the sensitivity is indeed dependent on the coherence time, i.e. the longer the coherent time is the 
better the sensitivity. In this section, we use the simple example of the radical pair with only isotropic hyperfine 
couplings to show that this is not the case in the present model. 



Mo+ = Tr 



a Tr 
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Figure 3: Decay of the entanglement in a radical pair reaction [Py-h'-[^ DMA-h'^^] under different types of quantum control. (N) 
without control; (Z) under Z control, (X) under X control. The curves are the estimated best lower bounds from the singlet 
fidelity, the symbols denote the values from numerical simulation. The magnetic field is B = 4.5 mT, and the control time is 
T = 0.5 ns. 



As we have described in the main text, for the Z control, we dynamically decouple the xx and yy hyperfinc inter- 
actions while keeping the magnetic-field dependent Zeeman interactions. Nevertheless, the magnetic field sensitivity 
is still much suppressed. This phenomenon can be understood as follows. The residual hyperfine couplings along the 
longitudinal direction (i.e. zz hyperfine couplings) only induce the transitions between the singlet state |S) and one 
specific triplet state |To), while these two eigenstates are degenerate and their energies are independent of the mag- 
netic field. In this case, the singlet-triplet interconversion is actually not influenced by the magnetic field, the effects 
of which are thus not detectable through the singlet yield. We will prove, in the following, that a chemical compass 
will loose its function, if one uses general dynamical decoupling protocols to promote the electron spin coherence. 

Assume that, at time to, the electron spins and the surrounding nuclear spins are in some state pitg) = pQ. The 
activation yield during a short time interval [to, to + t] is 



$(to,r)= r,{t)Ut)dt 

Jto 

with the singlet fidelity f{t) — {§\ps{t)\§) . We then write its first derivative with respect to the magnetic field as 

a$(to,to + r) 



(9) 



A(io,T) 



dB 



rc{t) at 



to 



dB 



which obviously determines the ultimate magnetic field sensitivity as A = -^("^'''i ''') summing up A(to,T) 
for all time intervals [to,to + T] = [mr, (m + l)r], m ~ 0,1,2,.... The singlet fidelity at time t G [to,to + t] is 
fs{t) = Tr (e~*^*^poe*^*^'Pal with At = t — to, and H is the Hamiltonian of the electron spins together with the 
nuclear spins as in Eq. (LI) [8|. By using a perturbation expansions for small At, we have 



= I - iAtH - IM!i/2 ^ o((At)3) 



(10) 



which enables us to express the singlet fidelity as follows. 



fs{t) = Tr 



po Vs+iAt[H,Vs] 



{Aty 



[[H,Vs],H] 



OiAt^) 



(11) 



where [A, B] = AB — BA. Using the properties of the singlet state that = 'Ps^ = 0, where ^ 

(2) 

Sz ), the first derivative of fs{t) can be written as 



(1) 



dfsjt) _ [Atf 
dB 2 



■Tr 



Po 



dH 

as 



HVs + VsH 



'dB 



0{{Atf) 



(12) 
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U Po = '-'^^ '^^^ easily verify that Tr [po (jfsHVs + ''^s-ff §§)] = 0. Thus, if a dynamical decoupling protocol is 

to protect the electron spin coherence during the reaction, i.e. keep the spin state close to the singlet state, we can con- 
clude that dfs{t)/dB ~ 0{{At)^), and A{to, r) will be of the fourth order in r, which is an order smaller than the one 
from other general states. We remark that we do not trivially assume that the system dynamics is frozen by protecting 
coherence, but the electron spin state does evolve even if it is kept closer to the singlet state under decoupling controls. 




Figure 4: Statistic of the visibility V for the initial radical pair state randomly chosen among the general product states (a) 
and the incoherent states (b), compared with the visibility for the singlet state Vs. 

Magnetic field sensitivity from random separable or incoherent states. — Here we compare the magnetic field sensitivity 
obtained from an initial singlet state (entangled) with the sensitivity obtainable from classically correlated states 
(separable) . A mixed quantum state is called separable if it can be written as a sum of product states 

P = ^Pk\(f>k)a{(l)k\<^\lpk)b{'lpk\ (13) 
k 

In general, a separable state can exhibit "coherence" , by which one means that some of its off-diagonal density matrix 
elements (with respect to the standard basis |t), arc non-zero. By definition, however, a separable state is not 
entangled. So there is an essential difference between entanglement and coherence. 

We introduce the optimal magnetic field sensitivity for the radical pair reaction [Py-h'^ DMA-h'^] on the set of 
separable states as 

Asep(S) = max |A(p,B)| (14) 

pGSep 

where A(p, B) denotes the magnetic-field sensitivity for a given initial state p. It can be proved that the optimal 
sensitivity AsepiB) is obtained by the product states |0)a ® \il')b- In our numerical calculations, we randomly choose 
5000 product states and calculate Asep{B). We also randomly choose the initial state from the set of incoherent 
states, the off-diagonal matrix elements of which are all zero, meaning that no coherence is present (and naturally 
no entanglement either). This is what we have done to compute the curve 'Sep' in Fig. Ll(b). We find, somehow 
surprising, that the optimal sensitivity from the incoherent states Ajnc is the same as Asep- In this sense, quantum 
coherence is not vital for the magnetic-field sensitivity. 

In the example of FADH'-O'", we characterize the angle dependence by using the quantity of visibility defined as 
follows 

max cf>i. - min $s 

^ = ^ ■ — T- (15) 

max <l>s + mm $s 

where $s is the singlet yield. We have randomly chosen a few hundreds of product states and incoherent states as 
the initial radical pair state, and calculated the corresponding visibility V (compared with the visibility Vs for the 
singlet state). It can be seen from Fig. |4]that a substantial part of separable (incoherent) states can account for an 
angular dependence that is as high as (or even higher) than for the singlet state. In this sense, the radical pair initial 
state of an avian compass need not be the singlet state. 
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Quantum entanglement and magnetic field sensitivity. — To quantify the amount of entanglement that exists in the 
active radical pairs during the reaction, similar to the singlet yield, we define the effective entanglement as the 
integral 



$B= / r,{t)E{t)dt (16) 
Its first derivative with respect to the magnetic field is the entanglement sensitivity 

A.= ^ (17) 

In Fig. [51 we plot A^; as a function of B. It can be seen that A^; changes conspicuously (kink) during the crossover 
between the regions of I and II. At the same time, the entanglement yield always increases with the magnetic field. 
This can be understood from the fact that strong magnetic fields will energetically suppress the relaxation (spin fiips) 
in the longitudinal direction. By this process, the state of the electron pairs changes towards a binary mixture of two 
entangled states (namely |§) and |To)), which is entangled for almost all values of the mixing parameter, resulting in 
a much longer lifetime of entanglement. 
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Figure 5: (Color online) Entanglement sensitivity A_b of a radical pair reaction [Py-^io DMA-ftj^^] as a function of the magnetic 
field B. The recombination rate constant is fc = 5.8 x 10*s~'^ 



To further illustrate the connection between quantum entanglement and the magnetic field sensitivity, we plot 
in Fig. [S] the time evolution of the entanglement and of the value of the accumulated magnetic-field sensitivity 
K{B,t) = d^{t)/dB, for different values of the magnetic field: B = 3mT, 3.5mT, 4mT, and 4.5mT. The hfetime of 
entanglement in the region of I is approximately Te = 4ns, while A(_B, t) needs about Tr = 10ns to reach its saturate 
value, see Fig.[6ja-b). We can also explicitly see the sudden increase of Tc when B crosses between the regions I (low 
magnetic field) and II (high magnetic field), from T^; = 4ns to about 7.3ns, see Fig. injc-d), which gives rise to the 
steps in Fig. L2 (b) Q. 

At this point, it is worth to emphasize that the concept of entanglement is different from the singlet fraction 
(fidelity), which was studied earlier. Generally speaking, a state can exhibit significant classical spin correlations 
without having any entanglement. In our specific example, at any time during the radical-pair reaction there will be 
a finite singlet fraction while entanglement, in contrast, will only exist for a much shorter time (as shown in Fig. 
In this sense, entanglement is a different, and generally more sensitive, signature than the singlet fraction. See also 
Fig. L2(b) in the main text. 

Reference model of a boson ic heat bath. — Let us assume that each of the unpaired electron spins is coupled with an 
independent bosonic heat bath at the same temperature. The dynamics of one central spin would thus be described 
by the following Lindblad type master equation 0, [l3| 

= ^i\H,,, p] + ^(2L,.pit _ pLlLk - lIl.p) (18) 



where He — rrihSz, with mt = — 7e_B, Li = ^/7S(t+ and L2 ~ ■\/7(l — s)cr_. The solution of the above master equation 
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Figure 6: (Color online) Time evolution of the (accumulated) magnetic-field sensitivity A{B,t) (rescaled) and the entanglement 
of the radical pair for (a) B = 3mT, (b) 3.5mT, (c) 4mT, and (d) 4.5mT in the radical pair reaction [Py-h'{^r)MA-h'{^]. The 
recombination rate constant is fc = 5.8 x 10* s^^ 



can be represented by a map p{t) ~ A^t[p(0)] which is explicitly expressed as follows 



Mt : It) (tl ^ at It) (tl + (1 - at) \i) {i\ 

^ (i-/3t)it)(ti+Ai;>ai 

It) ai ^ e-'2™,t m 
It) (tl ^ r,t It) (tl 



where at = (1 — s)e~^'>'*-|-s, /3t = se~'^'^^ + {l — s) and rjt = e~'^*. This map describes spin-exchange interactions with the 
environment with an effective rate 7 and an equilibrium parameter s that is related to the environment temperature 
T. The dependence of 7 and s on T and the magnetic field B is given in the following way: 7 ~ 2mi,Ko{2J\f + 1) and 
s — N'/(2Af + 1), where kq depends on the system-bath coupling strength on resonance, and the bosonic distribution 
function is Af = l/{e~ — 1) with the system energy scale eg = 2hmi, and the thermal energy scale = ki,T. Thus 
we have 



1 as 

1 57 
7^5 



e'T 

B tT 



1 



2 — e^T (e =T 



(19) 
(20) 



We are interested in the effects of low magnetic fields, for example B ~ \ mT, which corresponds to the thermal 
energy scale at temperature T ~ 2.69 mK that is quite low for biochemical systems. Thus we can naturally assume 



that ^ <^ \, from which it is easy to verify that 
smaller than 1 7 ^ | • 

The radical pair starts in the singlet state |§) = 



1 df 

7 dB 



\- — \ e E 



if T = IK then 



1 d-i 
7 dB 



is already four orders 



1 

72 



[WD ^ lit)) J and its state evolves as ps{t) — M. 



(1) 
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At time t, the density matrix is of the foUowing form 



where a = at{l - /3t), b ^ [atPt 
singlet fidehty fs{t) ^ Tr [p{t)Vs\ 



Ps{t) = 



/a \ 

c 

c 6 

V d y 



(21) 



+ (1 — at){l — Pt)] /2, d = (1 — at)l3t, and c = — ry^/2. Thus we can calculate the 

= — c as 



/,(t) = -[atA + (l-at)(l-A) 



The activation yield for the exponential re-encounter probability model is $ = fs{t)ke '^^dt, i.e. 



$ = 



k 



872 



k + 2"f (fc + 47)(fc + 27) 



.(1 



One can verify that under the general conditions we are interested in, the magnitude of the magnetic field sensitivity 
A would increase with the coupling strength scale kq, i.e. the fast thermalization is good in the present context. To 
achieve the optimal bound of A and illustrate the essential physics, we can assume that 7 is much larger than k (this 
is different from the real situation where 7 is smaller than fc), which leads to 



A~-(l-2s)— — 
B ex 



(22) 



The magnitude of A from the bosonic heat bath decreases as the temperature increases. Even at temperature as low 
as 1 K, it is already significantly smaller than the one from the nuclear spin environment, see Fig. [T] Therefore, we 
can conclude that the effects of low magnetic fields will indeed be washed out completely by the thermal fiuctuations. 
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Figure 7: (Color online) Magnetic field sensitivity A resulting from the nuclear spin environment (N) of the radical pair reaction 
[Py-/i^o DMA-/i^"5^]; and the optimal A (achieve when 7 3> fc) from the bosonic heat bath at temperature T = IK (see also Inset 
for an extended range of parameter) as a function of the magnetic field B. The recombination rate constant is fc = 5.8 x 10*s~^ 

i- 



By calculating d\A\/dB, we find that |A| will always grow as the magnetic field becomes stronger, as long as 
^ < ln(2 + \/3), which is obviously satisfied in the regions we are interested in. The change of the sign of d\A\/dB 
happens at |^ = ln(2 -I- -s/S) . At room temperature T = 300 K, this would correspond to the magnetic field B ~ 135T. 

The evolution of entanglement as obtained from Eq. ([21]) is E{t) = max{0, 2(|c| — (adY^^)}. In a similar way, 
one can obtain the lifetime of entanglement, see Fig. \8\ which is monotonically increasing with the magnetic field. 
This is another feature in marked contrast with the nuclear spin environment: there are no oscillatory kinks in the 
entanglement lifetime as the magnetic field increases. 
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Figure 8: (Color online) Lifetime of entanglement hqTe as a function of the magnetic field B. The bosonic thermal bath 
temperature is T = 1 K (red solid) and T = 1.5 K (blue dotted). 
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